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We review our recent theoretical works on the the quantum spin HaU effect. First we 
compare edge states in various 2D systems, and see whether they are robust or fragile 
against perturbations. Through the comparisons we see the robust nature of edge states in 
2D quantum spin Hall phases. We see how it is protected by the Z2 topological number, and 
reveal the nature of the Z2 topological number by studying the phase transition between 
the quantum spin Hall and insulator phases. We also review our theoretical proposal of the 
ultrathin bismuth film as a candidate to the 2D quantum spin Hall system. 

§1. Introduction 

Recently spin Hall effecP-l^ has been studied theoretically and experimentally. 
It has shed new light onto the physics of the spin-orbit coupling. In the spin Hall 
effect the spin-orbit coupling plays the role of the spin-dependent effective magnetic 
field, causing the Hall effect in a spin-dependent way. As a related subject, quantum 
spin Hall (QSH) phas^^'^S is recently proposed theoretically. The QSH phase 
is a novel topological phase, where the bulk is gapped and insulating while there 
are gapless states localized near the system boundaries. This phase can be found 
among nonmagnetic insulators, and opened up a renewed interest onto nonmagnetic 
insulators. This phase is a kind of topological order, and it is not evident compared 
with other type of orderings such as magnetism and superconductivity. This phase 
is not an ordered phase in the sense of Ginzburg-Landau theory. It is rather a 
topological order, which is encoded in the wavefunctions themselves, and it appears 
only at the boundaries. This topological order is hidden in the bulk, and appears as 
topologically protected gapless states at the sample boundaries and interfaces. The 
topological number, in the present case the Z2 topological number,'^ characterizes 
whether the system is in the QSH phase or not. In a sense, it plays the role of the 
"order parameter". This phase can be realized in 2D and in 3D®''^® the order 
guarantees the existence of gapless edge states for 2D and gapless surface states for 
3D. 

The QSH phase resembles the quantum Hall (QH) phase, while there are several 
important differences. The QSH phase requires an absence of magnetic field, while 
the QH phase requires a rather strong magnetic field. Moreover, the QH phase is 
realized usually in 2D, and it is not easily realized in 3D because the motion along 
the magnetic field usually does not become gapful. In contrast, in the QSH phase 
there is no such built-in direction and it can be easily realized in 3D as well as in 2D 
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without applying fields. 

These gapless edge states in 2D are peculiar in the sense that it is robust against 
nonmagnetic disorder and interactionP'^^^ This is in strong constrast with usual 
states localized at the boundary, which are sensitive to the boundary roughness and 
impurities and so forth. 

Since general readers may not be familiar with this issue, in the first half of this 
paper we give a general instructive review for the whole subject: "quantum spin Hall 
phase for pedestrians" . In the latter half we explain our recent research on this topic. 
The paper is organized as follows. In Section 2 we consider edge states of various 
2D systems and see whether they are robust or not. Surface states on the 3D QSH 
phases is also mentioned. In Section 3 we see how the quantum phase transition 
between the QSH and insulator phases occurs. Section 4 is devoted to explanations 
of an existence of gapless helical edge states based on the models obtained in Sec. 3. 
In Section 5 we theoretically propose that the bismuth ultrathin film will be a good 
candidate for the 2D QSH phase. In Section 6 we give concluding remarks. 

§2. Edge States in Various Systems — Fragile or Robust ? 

2.1. Edge States in Graphene 

Graphene has been studied intensively in recent years. One of the novel proper- 
ties of graphene is the edge states. It was theoretically proposed by Fujita et al.^ 
and was observed by STM experiments. This is a good starting point for studying 
the edge states in various systems. 

One simple way to see the edge states in graphene is the nearest-neighbor tight- 
binding model on a honeycomb lattice (Fig. HJa)), described as 

H = tY,4c,. (2-1) 

{id) 

We ignore other details of graphene, since they are inessential to the subsequent 
discussions. The primitive vectors are 

ai = ^(l,V3), a2 = ^(-1,^/3), (2-2) 

and the reciprocal lattice vectors are 

^ 27r / 1 \ _ 27r / 1 \ , , 

The Brillouin zone is shown as Fig. ^h). 

When we consider the tight-binding model on an infinite system, the eigenener- 
gies are 

E(k) = ±t e*''-^! + 6**^-^2 ^ I ^ (2-4) 

The gap between the valence and the conduction bands closes at two different points 
in the Brillouin zone, which are called K and K' points. Their wavenumbers are 

kK = ^(Gi - G2) = 0) , kK' = -kK. (2-5) 
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Fig. 1. (a) Honeycomb lattice and (b) Brillouin zone corresponding to the honeycomb lattice. 



Around these points the dispersion is hnear and forms Dirac fermions. To study 
edge states, we consider the system geometry with edges. For this purpose we make 
the system to have the ribbon geometry, having a finite width in one direction and 
being infinite in the direction perpendicular to it. The honeycomb lattice allows 
various types of edge shape, and we choose the zigzag edge and the armchair edge 
for example. The band structures for the two choices are shown in Fig. [2] (a) and (b), 
respectively. To understand the obtained band structures in this ribbon geometry, 
we relate the bulk band structure and that of the ribbon as follows. In the ribbon 
geometry, the translational symmetry in one direction (perpendicular to the ribbon) 
is lost and the wavenumber along this direction is no longer a good quantum number. 
Therefore the bulk band structure is projected along this direction. This almost 
corresponds to the band structure for the ribbon geometry, calculated in Fig. [2]^a) 
and (b). In the zigzag-edge case (Fig. [2^a)), however, the states located at = 0, 



27r 



< kxa < 



4tt 
3 



is outside the bulk-band projection. This is nothing but the edge 



states, because it is located in the bulk band gap, i.e. within the region where the 
bulk states are prohibited. 

We note that the above tight-binding model is with drastic simplification. In 
reality there are many factors which have been ignored here. We thus ask ourselves 
whether the above properties survive perturbations. For example, within the above 
tight-binding model, we can introduce various perturbations which do not support 
edge states. When we vary the boundary conditions, the edge states may vanish. 
For example, the armchair edge does not support edge states. Another perturbation 
which kills edge states around = is a staggered on-site potential, although it 
may not be realized in real systems. 

This fragility of the edge states arises because the edge states in graphene are not 
topologically protected, due to vanishing of the bulk gap. In the following sections we 
see various kinds of robust edge states. All of these robust edge states are associated 
with respective topological numbers. A more physical reason why the edge states in 
graphene are fragile is that no current of any kind is flowing along the edge. As we see 
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Fig. 2. Band structure of the nearest-neighbor tight-binding model of the the graphene in the 
ribbon geometry with (a) zigzag edges and (b) armchair edges. 



from Fig.[2]^a), the edge states form a flat band, which means that the vefocity afong 
the edge is zero. The modes are locaUzed and does not move along the edge. In fact, 
this flatness of the edge-state band is not a universal property; some perturbations 
give a dispersion to the otherwise flat band. For example, the next-nearest-neighbor 
hopping brings about a downward shift of the dispersion around kx = t^/cl, compared 
with Fig. [2] (a). Even in this case, the net current along the edge sums up to zero. 

2.2. Edge States in Quantum Hall Systems 

In this section we consider an example of robust edge states. We consider the in- 
teger QH systems. The QH systems are realized experimentally in a two-dimensional 
electron system in a strong magnetic field. In the QH systems the bulk is gapped 
while the edge has gapless edge states, which carries chiral current. The number of 
chiral edge states, u, is a topological quantity, which does not change under weak 
perturbation. In this case the robustness of the edge states comes from the topologi- 
cal number u. The bulk states bear a topological order, described by the topological 
number v called the Chern number. 

The edge states of the QH system look very different from those in the graphene 
discussed in the previous section. Nevertheless, by using a simple model we can relate 
them and discuss their differences. It is the model proposed by Haldane.^-' It is a 
tight-binding model on the honeycomb lattice, where the nearest-neighbor hopping 
is real, and the next-nearest neighbor hopping is complex. The Hamiltonian is given 
by 

Here Vij = sgn(di x d2)z = il, where di and d2 are unit vectors along the two bonds, 
which constitute the next-nearest neighbor hopping, represents a staggered on- 
site potential, and takes the values ±1 depending on the i-th sites being in the A 
or B sublattices, respectively. This simply means that the next nearest neighbor 
hopping, going around the hexagonal plaquette in the clockwise (counterclockwise) 
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way, obtains the phase e*"^ (e"**^). 

The quantization of axy for this model without impurities can be seen from the 
Kubo formula, and the resulting axy is rewritten as Uxy = ue^ /h^^'^^ where 



E (2-7) 
R7 — ^vr 



and 

_ _ ^«">(k) = (u„,| A |„„,) . ,2.8) 

The sum in Eq. (j2-7p is taken over the filled bands. The quantity has novel 
properties arising from topology. A naive application of the Stokes theorem casts 
Eq. (j2-7p into a contour integral along the border of the Brillouin zone (BZ), and 
the periodicity of the BZ results in = 0; however it is not true.'^ In some cases as 
the present one, the Bloch wavefunction cannot be expressed as a single continuous 
function over the whole BZ; the BZ should be divided into pieces, on each of which 
the Bloch wavefunction is continuous. At the borders between two "pieces" the 
Bloch wavefunctions differ by a U(l) phase. The number v is expressed in terms 
of this phase difference. As a result one can show that this quantity is quantized 
to be an integer, and it turns out to be the topological number called the Chern 
number.lIS'lHJ 

This Chern number represents the number of chiral gapless edge states going 
around the sample edge. Namely the existence of the gapless edge states is guar- 
anteed by the Chern number. This comes from the Laughlin argumentf^ we roll 
the two-dimensional system into an open cylinder by attaching two edges on the 
opposite sides, and let a flux ^ penetrate the hole. The flux <P is to be increased 
from to a flux quantum. Then the number of electrons carried from one edge of 
the cylinder to the other is equal to the Chern number. These carried electrons are 
on the gapless edge states. Thus we can establish the correspondence between the 
number of chiral edge states and the Chern number. 

To calculate the value of for the Haldane's model, it is not easy to calculate 
Eq. (|2-8p directly, and analytical calculation looks almost impossible. Instead by 
considering the change of by a change of some parameter, we can easily calculate 
u. First we notice that for = the system is time-reversal symmetric, which yields 
u = 0, because u changes sign under time-reversal. As long as the gap remains 
open, the Chern number u is quantized to be an integer and cannot change when a 
parameter in the Hamiltonian is continuously changed. In the Haldane's model, by 
changing the system parameters, the gap closes only at K and K' points. At the 
gap closing the spectrum of Eq. (j2-6p becomes linear in k near the gap-closing point, 
namely the massless Dirac fermion is formed at the gap closing. 

It is straightforward to show the following statement for the change of the Chern 
number at the gap closing. We consider a 2 x 2 Hamiltonian matrix H{M, k), which 
depends on a parameter M. Suppose the gap between the two bands closes at 
M = M(°) and k = k(°). Then we have F(M(o),k(°)) = £;o(mW, k(°))i, where i is 
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an identity matrix, and we can write 

HiM, k) = Eo{M, k)i + k)si, (2-9) 

i 

where Si are Pauh matrices. For notational brevity and convenience we write = M, 
ki = kx, k2 = ky. The coefficients Oj is expanded in M and k in the hnear order as 

a^= Yl ikj-kf^)aij. (2-10) 
i=o,i,2 

Then one can show that the change of the Chern number u across IS 

i^{M = m(°) + 6)- u{M = m(°) -S)= sgn(deta). (2-11) 

where a is the matrix with elements Uij and 6 is an infinitesimal positive number. 
Thus the Chern number changes by one at the gap closing. In the present case, the 
gap closes when M = =p3\/3t2 sin with k = ik/^- = (kx, k^/), respectively. For 
the respective cases, the matrix a is calculated by linearlizing the Hamiltonian in 
the vicinity of the gap closing as 

T\/3a/2 I (2-12) 

and sgn(deta) = ±1. Thus when we increase M across the value =p3\/3t2 sin the 
Chern number increases by ±1. From this we can easily elaborate the phase diagram 
as shown in Fig. [3j As we have seen, for analytical calculation it is much easier to 
calculate the change of the Chern number, than to calculate the Chern number itself. 




M 
-3/3 

Fig. 3. Phase diagram of the Haldane's honeycomb-lattice modeL f represents the Chern number 

What is the fundamental differnce between the QH phase and the grephene? The 
difference comes from the bulk gap. The graphene is gapless while the QH phase 
has a bulk gap. The topological number in the QH phase protects the existence of 
the gapless edge states. 
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2.3. Edge States in 2D Quantum Spin Hall Systems 

The QSH system is insulating in the bulk, while there are gapless edge states 
carrying a spin current. The simplest case of the QSH system is realized by super- 



{c^xy = —e^/h). The superposition of these states is the QSH state. The edge states 
consist of two states with opposite spins and velocities. These states are thus car- 
rying a spin current. To realize this state we need a spin-dependent magnetic field, 
which can be produced by the spin-orbit coupling. This example is just a superpo- 
sition of two QH systems with conserved spin . Nevertheless, in real systems, the 
spin-orbit coupling does not necessarily conserve s^, and the above two QH subsys- 
tems are mixed with each other. The next question is what happens when s^ is no 
longer a good quantum number. We consider this question in the following, and we 
see that the physics coming from topology survives partially. 

This phase is realized in the model proposed by Kane and Mele.l^® The model 



represents a staggered on-site potential, taking values ±1 depending on the sub- 
lattice index, djj is the unit vector along the nearest neighbor bond from i to j. 
For a special case, if A/? = and = 0, the model conserves s^ and it reduces to a 
superposition of two Haldane models: 



Hkm{Xr = 0, a. = 0) = i?i,idane('^ = -^/2) + ^Lldane('A = Vr/2). (2-14) 



For generic cases where s^ is not conserved, the key ingredient is the time-reversal 
symmetry, which gives rise to the Kramers degeneracy between k and — k. In the 
theory the wavenumbers satisfying k = — k (mod G) play an important role. Such 
momenta are called the time-reversal- invariant momenta (TRIM), and are expressed 



as k = where ri=(„^„2) = i(niGi +71262) with ni,n2 = 0, 1 in 2B^^^ The 



where is the time-reversal operator represented as = iayK with K being complex 
conjugation, n^.k is the periodic part of the m-th Bloch wavefunction lying below 
Ep, and N is the number of Kramers pairs below Ep. This matrix t/;(k) is unitary 
at any k, and is also antisymmetric at k = Tj. Then for each TRIM we define an 
index 5i as 






SzCj + iXn ^ c|(s X dij)zCj + K ^ iic\ci. (2-13) 





(2-16) 



for X-asymmetric systems,'^^ and 



N 




(2-17) 



m=l 
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for T-symmetric systems}^ where ^2m(^i) (= il) is the parity eigenvalue of the 
Kramers pairs at Fj. The index 5i takes the values ±1. The Z2 topological number 
I' is then defined as 

4 

(-i)^=n'^., (2-18) 

1=1 

where the product is taken over the TRIMs k = Fj. Hence (—1)'^ takes only two 
values ±1, which means there are only two distinct cases u = even and v = odd. 
For simplicity of notations we henceforth call these cases as = and 1^ = 1. If the 
resulting value is = 1 it is the QSH phase, and if it is = it is the ordinary 
insulator. The expression of the Z2 topological number is different between systems 
with and without X-symmetry. This shows a crucial role of the X-symmetry in the 
theory of the QSH systems. It is because the X symmetry is the only symmetry 
operation which relates between k and — k. 

To illustrate the physics of the Z2 topological number, we calculate the band 
structure for geometries with edges. Let us consider a ribbon geometry, which is 
finite in one direction and is infinite in the other direction. To see the difference 
between two phases, QSH and ordinary insulator (I), we take up the sets of the 
parameter values for the respective phases, employed in Ref. Hj). The result is shown 
in Figs, m and [5] for the QSH and I phases, respectively. As we see from Fig. HI there 
exist gapless edge states in the QSH phase, irrespective of the geometry. In contrast, 
for the insulator phase (Fig. ED there are no gapless edge states. In fact there are 
edge states, but they do not go across the gap. These edge states may or may not 
cross the Fermi energy. Even if they cross the Fermi energy, it is not an intrinsic 
property, and the crossing can disappear by perturbation. 




Fig. 4. Band structure of the Kane-Mele model in the QSH phase in the ribbon geometry with 
(a) zigzag edges and (b) armchair edges. The parameters are A„ = O.lt, Xso = 0.06t and 
\r = 0.05t. 

The bulk band-structure is gapped for both the QSH and I phases and looks 
similar. Namely, the topological order is not evident in the bulk band structure. 
The existence of the robust edge states is encoded in the bulk wavefunctions. There 
is a bulk-edge correspondence: 

1. Z2 topological number is = 1 (i^ = 0). 
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Fig. 5. Band structure of the Kane-Mele model in the I phase in the ribbon geometry with (a) zigzag 
edges and (b) armchair edges. The parameters are \v = QAt, \so — 0.06t and \r = O.OSt. 



2. there are odd (even) number of Kramers pairs of gapless edge states. 
This correspondence can be shown by the argument similar to the weU-known Laugh- 
hn's gedanken experiment.'^''^ Suppose we consider the system on a ribbon, with 
two opposite ends attached. In one direction the system is periodic while in the other 
directions there are edges. When we increase the flux penetrating the hole from zero 
to half of the flux quantum, the change of a certain physical quantity (time-reversal 
polarization) is zero for v = ^ while it is unity for v = 1."^ This means that for 
V = \ there is a Kramers pair of gapless edge states while for = there are no 
gapless edge states. 

2.4. Surface States in 3D Quantum Spin Hall Systems 

The analogous phase is also possible in 3D. In this case this phase is an insulator 
in the bulk and supports gapless surface states carrying spin currents. In this case 
as well, there is a correspondence between the bulk and the surface. The topology 
of Fermi "curve" of the surface states is related with the Z2 topological numbers for 
the bulk. 

One can see this by the 3D tight-binding model introduced by Fu et al^ on a 
diamond lattice. This model exhibits a transition between QSH and I phases. The 
model is written as 

H = tY, c\cj + i(8Aso/a') • (4 (2-19) 

fe-> m) 

Here a is the size of the cubic unit cell, t is the nearest-neighbor hopping, and 
s = {sx,Sy,Sz) are the Pauli matrices. The term with Aso is a spin-dependent 
hopping to the next nearest neighbor sites, representing the spin-orbit coupling. 
The vectors d^- and d^^ denote those for the two nearest neighbor bonds involved in 
the next-nearest-neighbor hopping. 

In 3D, the TRIMs are rj=(„j„2n3) = \{niGi + n2G2 + n^G^) with ni,n2,n3 = 
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0, 1. There are four Z2 topological numbers 1^0,1^1,1^2, u^^^ given by 

8 

(-ir=n<5., (-1)'^'== n ^^=in^n,n,)■ (2-20) 

1=1 nfe=l;raj_^fc=0,l 

Each phase is expressed as lyo; {1^11^21^3) , which distinguishes 16 phases. Because 
among Vi, uq is the only topological number which is robust against disorder, the 
phases are mainly classified by vq. When uq is odd the phase is called the strong 
topological insulator (STI), and when i/q is even it is called the weak topological in- 
sulator (WTI). The STI and WTI correspond to the QSH and I phases, respectively. 
The other indices vi, f 2 , and v^, are used to distinguish various phases in the STI or 
WTI phases, and each phase can be associated with a mod 2 reciprocal lattice vector 
G'i/ii/2V3 = '^ibi + z^2b2 + ^'3b3, as was proposed in Ref. [6]). These topological numbers 
in 3D determine the topology of the surface states for arbitrary crystal surfaces.'^ 
We note that among the four Z2 topological numbers in 3D, only is robust against 
nonmagnetic impurities, while the others {v^ (k = 1, 2, 3)) are meaningful only for a 
relatively clean sample.^^ 

§3. Quantum Phase Transition with a Change of Topological Number 

As we have seen in the Haldane's model on the honeycomb lattice, it is hard 
to calculate the topological number itself and to capture its physical meaning in an 
intuitive way. The topological number involves an integral over the whole Brillouin 
zone, both for the QSH systems and the QH systems. 

On the other hand, the "change" of the topological number is more accessible. 
This is because the change occurs locally in k space. As we change an external 
parameter, the system may undergo a phase transition. It accompanies a closing 
of a bulk gap at a certain k, because it is the only way to change the topological 
number. When an external parameter is changed, in some cases the phase transition 
occurs, while in other cases it does not, because of the level repulsion. 

In the following we consider "generic" gap closing by tuning a single parameter, 
which we call m. We exclude the cases where the gap closing is achieved by tuning 
more than one parameters. In such cases, the phase transition might be circum- 
vented by some perturbation. In general, energy levels repel each other, thereby 
the valence and the conduction bands do not touch if the number of tuned param- 
eters is not large enough. The number of tuned parameters to achieve degeneracy, 
called the codimension, is sensitive to the symmetry and the dimension of the system 
considered. 

We henceforth consider only clean systems without any impurities or disorder. 
The time-reversal symmetry is assumed. We also assume that the Hamiltonian is 
generic, and we exclude the Hamiltonians which require fine tuning of parameters. 
In other words, we exclude the cases which are vanishingly improbable as a real 
material. 
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3.1. 2D Quantum Hall Systems 

For the 2D QH systems, it is simple, and in fact we have already seen the 
physics of gap closing in the Haldane's model. We can argue gap closing in generic 
systems, and it becomes similar to the one already described in Section 2. The Chern 
number is the total flux of Bz(k) inside the Brillouin zone. Therefore, if we roll the 
Brillouin zone into a torus, the Chern number is nothing but a total monopole charge 
inside the torus, where a monopole and an antimonopole have charges +1 and — 1, 
respectively. Each band is associated with the respective Chern number. At the gap 
closing, the Chern number at the lower band changes by ±1 and that of the upper 
band changes by =Fl- This means that a monopole (or an antimonopole) goes out of 
the Brillouin-zone torus of the upper band, and into the Brillouin-zone torus of the 
lower band. 

Let m denote the parameter which drives the phase transition. In the Hal- 
dane's model the on-site staggered potential M plays the role of the parameter m. 
The feature of the phase transition is further clarified by considering a hyperspace 
(m, kx,ky) = {ko, ki, k2), and characterizing the gap closing in terms of the gauge 
field in m-k.-space^^'^^ Suppose the gap closes at an isolated point k = (m, k). 
Then the involved bands, which we call a-th and /?-th bands, have monopoles in 
m-k space, with their monopole charges are opposite in sign.'I^'l^ More precisely, 
the gauge field (k) and the corresponding field strength Bq (k) for the a-th band 
are defined as 

A«(k) = -i(V'«(k)|VklV'«(k)), (3-1) 
B,(k) = X A„(k). (3-2) 
The corresponding monopole density is defined as 

/5„(k) = i-V^-B<,(k) (3-3) 

Except for the point where a-th band touches with other bands, the monopole den- 
sity /9a (k) vanishes identically. At the k point where the a-th band touches with 
another band (/3-th band), the wavefunction cannot be written as a single analytic 
function around this point, and the wavefunction is to be written with more than 
one "patches" which are related by gauge transformatioUj^SJ as is similar to the vec- 
tor potential around the Dirac monopole in electromagnetism.'^ This leads to a 
(5-function singularity of p(k) at the band touching; /9a(k) ~ — /9^(k) ~ ±d(k = ko). 
As a result the monopole density is written in general as p(k) = gz5(k — k;), where 
qi is an integer representing a monopole charge. The monopole charge is conserved 
under a continuous change of the Hamiltonian. The monopoles indicate the value 
of (m, k) where the gap closes. At the monopole, the Chern number for the band 
considered changes by ±1. 

As an example, in the Haldane's model on the honeycomb lattice, the gap closes 
either at K or at K'. When the gap closing occurs either at K or at K' , the Chern 
number changes by one. Meanwhile, when M is changed while (j) = is kept, the 
gap closes simultaneously at K and K', and the changes of the Chern number at K 
and K' cancel each other; the Chern number does not change as a result. 
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3.2. Quantum Spin Hall Phase and Universal Phase Diagram 

Because the QSH phase is roughly analogous to a superposition of two QH 
systems, the phase transition can be studied similarly. The Z2 topological number 
is preserved as long as the gap remains open. Suppose the system goes from the 
insulator to the QSH phase by changing some parameter of the system. Then the 
gap should close somewhere in between. Whether or not gap closes as the parameter 
is varied reflects the topological properties of the system. We investigate the criterion 
for the occurence of the phase transition in Refs. I27p . l28|) .l29 p . as we explain in the 
following. Through this study we can show that the gap-closing physics is equivalent 
to the physics of the Z2 topological number. 



(a) 3D 



gapless 



symmetric 




(b)2D 




symmetric 



I \ -QSH 



Fig. 6. Universal phase diagram between the QSH and the insulator phases in (a) 3D and (b) 2D. 
m is a parameter driving the phase transition, and 5 is a parameter describing the degree of 
inversion-symmetry-breaking. 

To study the phase transition in 2D and in 3D, we consider a Hamiltonian matrix 



(3-4) 



We assume that the system is a band insulator, and the Fermi energy Ep lies within 
the gap. The time-reversal-symmetry gives, 



H{k) 



SyH 



\-'k)Sy, 



which is rewritten as /ii-|(k) = hj, {—]<.), /i||(k) 



-k), /liT(k) 



(3-5) 
-k). 



The Kramers theorem guarantees that the band structure of such time-reversal- 
symmetric spin-1/2 system is symmetric with respect to k <-> — k. While the dimen- 
sion of the Hamiltonian is arbitrary, it will be taken as the number of states involved 
in the gap closing. 

The feature of the phase transition is different whether the system considered 
is (i) T-symmetric or (ii) X-asymmetric.'^''^ It is because the degeneracy for each 
state is different for the two cases. According to the Kramers theorem, the time- 
reversal-symmetry says ena(k) = £nai—^), where ena(k) is the energy of the n-th 
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Fig. 7. Generic gap-closing for (a-1) 2D inversion-asymmetric, (a-2) 2D inversion-symmetric, (b-1) 
3D inversion-asymmetric and (b-2) 3D inversion-symmetric cases. In the cases (a-2) (b-2) all the 
states are doubly degenerate by Kramers theorem. In (a-l)(a-2) and (b-2), the gap closing and 
concomitant phase transition occurs only at a single value of m: m = mo. Meanwhile in (b-1), 
by increasing m, the gap closes at m = mi, and the bulk remains gapless in mi < m < m2. 
The gap opens again at m = m2. Although in reality the k space is two-dimensional in (a- 
l)(a-2), and three-dimensional in (b-l)(b-2), it is drawn as one- dimensional in (a-l)(a-2), and 
two-dimensional in (b-1) (b-2) for clarity of illustration. 

band with pseudospin a, and a is the pseudospin opposite to a. If in addition, 
the system is X-symmetric (i), all the states are doubly degenerate, because the 
T-symmetry imposes ena(k) = k), leading to ena(k) = e„Q(k). If (ii) T- 

symmetry is broken, double degeneracy occurs only at points k = Fj, where Fj is 
one of the TRIM. By analyzing the respective cases (i)(ii) in 2D and 3D, we obtain 
a universal phase diagram shown in Fig: [6l Here m is a parameter driving the phase 
transition, and a parameter 6 describes the degree of inversion-symmetry-breaking. 
The derivation of this universal phase diagram is generic and based on topological 
arguments. Hence it does not depend on the details of the system. The parameter m 
can be any parameter, and in the CdTe/HgTe/CdTe quantum well the well thickness 
d plays the role of m here. 

3.2.1. Inversion asymmetric systems 

In T-asymmetric systems, when k 7^ Fj, each band is non-degenerate. At the 
gap-closing point, one valence band and one conduction band become degenerate. 
In this case a 2 x 2 Hamiltonian matrix is sufficient for our purpose; 

«=(;. I). (3-6) 

where o, b are real functions of k and m, and c is a complex function of k and 
m. A necessary condition for the two eigenvalues to be identical consists of three 
conditions a = b, Rec = and Imc = 0, i.e. the codimension is three .'^^ '13311 To put 
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it in a different way, the 2x2 Hamiltonian i7(m, k) is expanded as 

3 

if(m, k) = ao(m, k) + aj(m, k)cji. (3'7) 

i=l 

The gap closes when the two eigenvalues are identical, i.e. when the three conditions 
aj(m, k)=0 (i = 1,2,3) are satisfied. This means that the codimension is three. 

In 2D, the codimension three is equal to the number of parameters involved, 
that is, kx, ky and m. Thus the gap can close at some k when the parameter m 
is tuned to a critical value. Near the gap-closing point k = ko(7^ Fj), the system's 
Hamiltonian corresponds to massive Dirac fermion, and can be expressed as 

H = {m- mo)az + (k^ - kox)crx + {ky - koy)ay (3-8) 

after unitary and scale transformations. The time-reversal-symmetry requires that 
the gap closes simultaneously at ko and — ko as shown in Fig. [7] (a-1), and that the 
masses of the Dirac fermions at k = itko have opposite signs. In the Kane-Mele 
model for the QSH phas^® the gap closes simultaneously at the K, K' points, 
corresponding to the present case. 

In 3D, as is different from 2D, the gap closing at k = itko 7^ Fj cannot lead 
to phase transion. This is because the codimension three is less than the number 
of parameters {m,kx,ky,kz)- The three gap-closing conditions determine a curve 
in the four-dimensional space {m,kx,ky,kz). When m is changed continuously the 
gap-closing k point moves in the k space, and the system remains gapless. This 
curve forms a loop C in m-k space, and the gap opens when m is changed across 
the extremum of the loop. The loop C is a trajectory of the gap-closing points i.e. 
monopoles in the k-space as we change m. It forms a closed loop in the m-k space, 
because of the conservation of the monopole charge. The time-reversal-symmetry 
requires 

B„(k) = -Bs(-k), /j„(k) = ps(-k). (3-9) 

The monopoles are symmetric with respect to the origin. Therefore, the generic 
form of the loop C is as shown in Fig. [Sj Thus the loop C occupies a finite region 
in the value of m, and it follows that in X-asymmetric 3D systems, gapless phase 
emerges, '251' which is nonexistent in 2D. 

So far we discussed the gap closing at k / Fj. To complete the discussion, we 
show that the gap does not close at k = Fj, when the system is X-asymmetric. At 
k = Fj, the band is doubly degenerate, and the codimension is five,!^''^^ exceeding 
the number of tunable parameters which is one (that is, ra). Thus, generic gap- 
closing cannot occur at k = Fj. One can see this as follows. Because both the valence 
and the conduction bands are doubly degenerate, we consider 4x4 Hamiltonian 
matrix with time-reversal-symmetry. From (j3-5p we get 



i=l 



(3-10) 
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Fig. 8. Trajectory of the gap-closing points for (a) inversion- (I-)asymmetric and (b) symmetric 
systems. For (b) T-symmetric systems, the gap-closing point is located at k = Fi, and isolated 
in the m-k space. Only at rn = mo the system is gapless. For (a) T-asymmotric systems, 
the gap-closing points are created in monopole-antimonopole pairs at m = mi, and move in 
k-space as m changes. Solid and broken curves denote the trajectories of the monopoles and 
antimonopoles, respectively. The system opens a gap when these gapless points annihilate in 
pairs at m = m2. 



where Oj's and Eq are real, and i^i = 1 (g) t^, r-2 = Gz ® Ty, 73 = 1 (g) r^, = 
(Ty <Si Ty, and Is = ax <Si Ty, and crj, Tj are the Pauli matrices. The eigenenergies 

are Eq it \jYlt=i ^f- "^^^^ condition for the gap-closing between the two (doubly- 
dcgcncratc) bands consists of five equations = for i = 1, • • • ,5, wliich arc not 
satisfied by tuning only one parameter m. (Here the wavenumber k is fixed to be 
Fj.) Thus the gap does not close at k = Fj by tuning a single parameter m. 

3.2.2. X-symmetric systems 

In T-symmetric systems, the energies are doubly degenerate for every k by the 

Kramers theorem. The gap closes between the two doubly-degenerate bands, and 
we set the Hamiltonian matrix H(k) to be 4x4. The T-symmetry is imposed as 

H{-k) = PH{k)p-\ u(-k) = Pu(k), (3-11) 

where P is a unitary matrix independent of k which commutes with the spin matrices, 
and u(k) is the periodic part of the Bloch wavcfunction: ip]i{r) = u(k)e^^"^ . Because 
= 1, the eigenvalues of P are ±1. By a unitary transformation which diagonalizes 
P, one can rewrite 

P=(^^^ P^ ) ' = n = diag(r/„, % = ±1, Vb = ±1 (3-12) 

without losing generality, rja and 77^ represent the parity eigenvalues of the wavefunc- 
tions. One of them corresponds to the valence band, and the other to the conduction 
band. 

While there are four combinations for r^a = ±1 and rji, = ±1, the overall sign 
for [rja, rib) is arbitrarily changed by gauge transformation. Thus the only distinct 
case are (i) r]a = rjb and (ii) rja = —rjb- The case (i) rja = rjb = ±1 means that the 
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wavefunctions (orbitals) a, b have the same parity, e.g. two s-hke orbitals or two 
p-like orbitals. When r]a = rjb = il; the Hamiltonian becomes 

5 

H{k) = Eo{k) + «i(k)ri, (3-13) 

i=l 

where Oj's and Eq are real even functions of k. On the other hand, when (ii) 
r]a = —rji) = ±1, where the two constituent wavefunctions have different parity, 
the Hamiltonian is 

4 

H{k) = Eo{k) + a,{k)n + 5^6»(k)i;' (3-14) 

1=1 

where -Eo(k) and a5(k) are even functions of k, and 6^*''(k) are odd functions of k. 
The matrices r{ = az ^ Tx, r2 = I (E) Ty, = ax(^ t^, r^ = ay^ Tx, and = 1 (gi 
form the Chfford algebra. Therefore, in both cases, r]a = r]h and rja = —rjh, the 
codimension is five, which exceeds the number of tunable parameters (m, kx, ky, 
kz). Therefore, gap closing does not occur in general. 

However, this counting holds true only when k is at a generic point with k 7^ Fj, 
At k = Fj, the codimension (number of parameters to achieve degeneracy) remains 5 
for r]a = r]b (Eq. (j3-13p ). while it becomes 1 for r]a = — % (Eq. ()3-14p ). This is because 
the odd functions 5^*-* (k) vanish identically at k = F j , and one has only to tune 05 (k) 
to be zero. The number of parameters (m) is one, because the wavenumber is fixed as 
k = Fj. Thus, the gap closes by fine-tuning a single parameter, only when ija = —r]b- 

In the following we derive an effective Hamiltonian describing the low-energy 
physics of the QSH-I phase transition. As a gap-closing point, we take k = as 
an example, and write down the Hamiltonian explicitly. Extension to other k = Fj 
points is straightforward. The Hamiltonian is expanded in terms of k as 



H{m, k) ~ So + + (/3« • k) (3-15) 



i=l 



where Eq and m are constants, and /J^*-* (i = I,-- - ,4) are two-dimensional real 
constant vectors. The critical value of m is set as zero. After unitary transformations, 
the Hamiltonian finally becomes block-diagonal. 



H{m, k)=Eo + 



/ m ^ 

— m 

m — z+ 
y — z_ —m j 



(3-16) 



where z± = h\kx -|- h-^ky it ih^ky with real constants 61, 62 and 63. If the system has 
fourfold rotational symmetry in the xy plane for example, one has h\ = 62 and 63 = 0, 
leading to z± oc kxUky. Thus we have shown that a generic Hamiltonian with time- 
reversal- and 2- symmetries decouples into a pair of Hamiltonians describing two- 
component Dirac fermions, with opposite signs of the mass terms. Such decoupling 
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is nontrivial. This Hamiltonian is identical with the one suggested for the HgTe 
quantum well in Ref. 24). The eigenenergies are E = Eoi: \J w? + and the gap 
closes at k = when the parameter m is tuned to zero. This kind of effective model 
can be used for studying disorder effects in the QSH phase's® 

3.2.3. Zi Topological Number 

We have discussed generic types of gap closing in time-reversal invariant systems, 
achieved by tuning a single parameter. There are two types of gap closing: (a) 
simultaneous gap closing at k = ibko(7^ Tj) occurs in systems without T-symmetry, 
and (b) gap closing between two Kramers-degenerate bands (i.e. four bands) at 
k = Fj occurs in systems with X-symmetry (see Fig. [7|). Thus the gap-closing by 
tuning a single parameter occurs only in limited cases, due to level repulsion. We note 
that we have not made any assumption on the topological number, in deriving this 
result. Nevertheless, we can show that the both cases (a) and (b) involve a change of 
the Z2 topolgical number, and are accompanied by a quantum phase transition. The 
Kane-Mele model on the honeycomb lattic^ belongs to class (a) while the HgTe 
quantum- well modeP^ belongs to class (b). 

In 2D T-symmetric systems, the gap closing at the QSH-I transition occurs 
at TRIM k = Fj. This is accompanied by an exchange of the parity eigenvalues 
between the valence and the conduction bands. It corresponds to the expression 
of the Z2 topological number as a product of the parity eigenvalues over all the 
TRIMs k = Fj over the occupied states^*^' (Eq. (pTT]) ). On the other hand, for 
X-asymmetric 2D systems, the gap closes at ibko(7^ Fj) by tuning m. Because the 
Z2 topological number should change at the gap closing, the Z2 topological number 
should be expressed as an integral over the k space. This is the Pfaffian expression 
of the Z2 topological number (Eq. (j2T6p ).l^'l^ In retrospect, the gap closing in the 
X-symmetric systems occur only at Fj between the bands with the opposite parities. 
This is because otherwise the codimension is five, exceeding the number of tunable 
parameters. On the other hand, if X-symmetry is absent, and the level repulsion 
is less stringent, and the gap can close at some k other than the TRIM Fj. This 
difference is reflected in the definition of the Z2 topological number. 

3.3. Example: 3D Fu-Kane-Mele model 

The 3D Fu-Kane-Mele modeP is an ideal model for studying QSH-I phase tran- 
sitions in 3D. It is a 4-band tight-binding model on a diamond lattice and is X- and 
time-reversal-symmetric. It means that every eigenstate is doubly degenerate by the 
Kramers theorem. The doubly-degenerate conduction band and the valence band 
touch at the three X points, X'^ = (27r/a)r (r = x,y,z). To describe the phases 
having a bulk gap, one sets the nearest-neighbor hoppings for four bond directions 
to be different; ti {i = 1,2,3, 4).!^ The system then opens a gap. When we set the 
hopping to be tj = i -|- 6ti and St^ = = dt^, the phase diagram is as shown in 
Fig. El[a) as a function of 5ti and 5t2, obtained in Ref. E]). At the phase boundaries 
the bulk gap vanishes. 

To verify the universal phase diagram in 3D, we introduce the X-symmetry- 
breaking term, which does not exist in the original Fu-Kane-Mele model. The sim- 
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(a) (b) 




5fi 6^^= 5^1+ 5^2 



Fig. 9. Phase diagrams for the Fu-Kane-Mele model with Sts = 0, 5t4 — 0. ti and t2 are the bonds 
along the 111 and 111 directions. We put Aso = O.lt. The axes are in the unit of t. (a) The 
phase diagram in St\-5t2 plane.'^ A„ is set as zero, (b) The phase diagram in the St+-\v plane. 
Here 5t+ = 5t\ + 5t2, while 5t^ — 5ti — 5t2 = O.lt is fixed. The arrows in (a) and (b) refer to 
the same variation of parameters. 



plest way to break T-symmetry is to introduce an alternating on-site energy A„ into 
the system, as was done the 2D Kane-Mele model on the honeycomb lattice.l^ 

We then calculate how the WTI-STI phase transition is modified by the A^, 
term. As we see from the phase diagram (Fig. O^a)) , we regard 5t+ = 5ti + 5t2 as a 
parameter m driving the phase transition, while 5ti — 5t2 is fixed to be 5ti — 5t2 = O.lt 
as an example. This corresponds to the arrow in Fig. [D^a). The phase diagram in the 
(5t+-A„ plane is calculated as shown in Fig. [9l^b). When the X-symmetry is broken 
(At, 7^ 0), the gapless region appears in the phase diagram, in accordance with our 
universal phase diagram.!^ 

The trajectory ("string") of the gapless points in k space also agrees with our 
theory. As the parameter (5t+ is changed along the arrow in Fig. [TUlfa). the gapless 
points move in k space as in Fig. [TOTb).^^- The overall feature of the trajectory, 
i.e. its pair creation and annihilation with changing partners, perfectly agrees with 
our theory. The change of the Z2 topological numbers is also consistent with our 
theoryP 

§4. Helical edge states 

From the effective model thus obtained, we can calculate the helical edge states 
appearing at the boundary between the QSH and the insulator phase. Such boundary 
is described by setting the mass parameter m{x) to be dependent on space; m(ibc>o) = 
ibmo, i.e., 

/ mo : X > , . 

m = < „ 4-1 

[ -mo : X < 0. ^ ^ 

The detail of the crossover between mo and —mo is unimportant and is left unspeci- 
fied. For 2D T- asymmetric systems (Fig. [7)^ a- 1)), one can consider the Dirac fermions 
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Fig. 10. Trajectory of the gapless points in k space. As we change 5t+ with A^, fixed as shown in 
the arrow in (a), the monopoles and antimonopoles travel in the k space as shown in (b) by the 
solid and broken curves, respectively. The wavenumber k is shown in the unit of (27r/a). 



at k = ibko separately. Masses of these Dirac fermions change sign at m = 0; hence 
they yield the edge states localized at the boundary, as explained in Ref. [38j) . Be- 
cause the Dirac fermions at k = itko are related by time-reversal-reversal symmetry, 
the two edge states form a Kramers pair and carry a spin current. 

For 2D X-symmetric systems (Fig. [7]^a-2)), we follow the discussion in Refs. 1551). 
39|) to show that such a boundary between phases with different Z2 topological 
numbers has a Kramers pair of edge states. By replacing kx by —idx in Eq. (I3-I6p . 
we consider 



H{ky) = Eo + bid, 

+ 

V 



f -i 
-i 

V 



\ 

i 

1 / 



m (63 - ib2)ky 

(63 + ib2)ky -m 



m -{b-i + ib2)ky 

-{bs-ib2)ky -m 



(4-2) 



To calculate the eigenstates it is convenient to perform unitary transformation as 



H'{ky) = Q^H{k)Q = Eq + 



( b2ky m - bidx 
m + bidx -b2ky 



where 



Q 



T2 



\ 



-ibskyx/bi 



—b2ky m — bid: 
m + bidx b2ky ) 



( 1 1 



(4-3) 
(4-4) 



-1 1 j 
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The eigenvalue problem reads as H' {ky)uky{x) = E{ky)uky{x). The term Eq is ab- 
sorbed by shifting the energy. Because ()4-3p is block-diagonal, we first solve the eigen- 
value problem for the first two components of u^y- By putting u^y = (ui, n2, 0, 0)*, 
we get 



{E - b2ky)ui 
{E + b2ky)u2 



DU2, 



(4-5) 
(4-6) 



where D 
respectively: 



m 



= rn + bi-^,. They yield eigenequations for ui and U2, 



- dx 



DD^ui 
D^Du2 



{E^ 
{E^ 



blkl)ui, 
blkl)u2. 



(4-7) 
(4-8) 



Because (|4-7p is invariant under E — > —E, the resulting spectrum seems to be 
symmetric with respect to = 0; —E. However, it is not true, because 

in some cases the ui solutions to (j4-7|) has no corresponding solution for U2- If 



-bok 



J/' 



E 

solved for ui. Thus the solutions which are not symmetric with respect to E 



^ cannot be solved for U2. Similarly, li E = b2ky, (j4-6p cannot be 






are as follows. For ui{^ 0) which satisfies D^ui = 0, we get E 
from Eqs. (j4-5p and (j4-6p . whereas there is no solution with E = 



■ b2ky and U2 
b2ky. In the same 

token, for U2 which satisfies Du2 = 0, we get E = —b2ky from (j4-5p . whereas there is 
no solution with E = 62 fc^. Hence the spectral asymmetry is related to the kernels 
for D and . For example, for 61 > and mo > 0, the solution at the boundary 
dm), with D^ui = gives 



ui oc exp 



m{s)ds 



(4-9) 



and E = b2ky, while Du2 



has no normalizable solution. Thus the energy disper- 
b2ky, which crosses the Fermi energy E' ~ 0. 
This state is gap less, localized near x = 0. We have thus far solved the eigenequation 



sion in ky direction has a branch E 



Fig. 11. Schematic dispersion curves for the model (l4-2p . 
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for the first two components. The lower two components of the wavefunction u is ob- 
tained from above by time-reversal operation. Therefore, the above-mentioned edge 
state with E = 62^?/ has a Kramers partner with E = —b2ky. The whole dispersion 
is shown in Fig. [TTl Thus we have shown that the Kramers pair of edge states exists 
at the boundary between the QSH and I phases. They cross at ky = 0, as follows 
from the Kramers theorem. 

§5. Bismuth Ultrathin Films 

The QSH phase requires no magnetic field. This means that some materials 
might realize the QSH by themselves without applying any field. The only necessary 
conditions for the QSH systems are as follows. 

1. nonmagnetic insulator 

2. the Z2 topological number is odd = 1) 

The latter condition means that the spin-orbit coupling should be strong enough, 
which requires relatively heavier elements. In the absence of the spin-orbit coupling 
the Z2 topological number v is zero (i.e. even). When the spin-orbit coupling is 
made stronger, some systems can change its Z2 topological number from = to 
u = 1. At the phase transition, the gap closes. 

In this sense, the gap should be opened by the spin-orbit coupling. This is an 
ambiguous statement. In Ref. l26p we clarified that systems with large susceptibility is 
a good starting point for the materials search. Among materials with heavy elements, 
we pick up bismuth as a candidate. Bismuth is known as a strong diamagnet, due 
to interband matrix elements by the spin-orbit coupling. In this sense the gap is 
originated from the spin-orbit coupling, and is a good candidate for the 2D QSH. 
Bismuth itself is a nonmagnetic semimetal, not an insulator. We have to open a 
gap by some means to make it the QSH phase. One idea is to make it into thin 
film. Indeed, the recent experiments and first-principle calculations show that the 
(111) 1-bilayer bismuth has indeed a gap.l^ In Ref. I26p we considered the (111) 1- 
bilayer bismuth ultrathin film from the simple tight-binding model, and theoretically 
proposed that it is the QSH phase. We also calculate the parity eigenvalues for the 
tight-binding model, and confirmed the result. 

The other way to open a gap is to make an alloy with Sb. This leads us to the 
3D QSH, as has been proposed in Ref. [6]). In the ARPES measurement ,^23 the Fermi 
"surface" for the surface states has been observed, and it crosses the Fermi energy 
odd times, from the F point to the M point. This shows that the surface state of 
Bio.gSbo.i is the topological one for the QSH system. 

§6. Concluding Remarks 

With simple examples we have seen various kinds of edge states. In graphene, 
the existence of edge states is sensitive to boundary conditions, while in the QH 
and the QSH systems, the gapless edge states exist irrespective of the boundary 
conditions. This comes from the nontrivial topological number carried by the bulk 
states, defined only for insulators. On the other hand, the graphene, being a zero-gap 
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semiconductor, cannot have such topological numbers, which means that the edge 
states are not robust against perturbations. 

The quantum phase transition between the QSH and insulator phases is stud- 
ied. We consider generic time-reversal-invariant system with a gap, and study the 
condition when the bulk gap closes by tuning a single external parameter. Due to 
level repulsion, the gap does not always close by tuning a single parameter; instead 
in many cases, fine tuning of more than one parameters is needed to close the gap. 
In the T-symmetric systems, the gap closes only at the TRIM k = Fj, between the 
valence and conduction bands with opposite parities. In the X-asymmetric systems, 
on the other hand, the phase transition is different between 2D and 3D. In 2D the 
gap closes simultaneously at k = itko / Fj. In 3D there appears a gapless region 
in the phase diagram between the QSH and the insulator phases. The gap closing 
points are monopoles and antimonopoles, and they are created and annihilated in 
pairs, when the system transits from the gapless phase into the phases with a bulk 
gap (i.e. QSH or insulator phases). 

It is interesting to note that in each case with robust gapless edge states, there 
is an associated current. In the 2D QH phase it is the charge current which the edge 
states carry. In the 2D QSH phase it is the spin current. There are also other classics 
of systems which show this kind of stable edge states. One is the superconductor 
without time-reversal symmetry, for example with gap function + ipy ■ In this case 
the current of Majorana fermion is carried by the edge states, and the edge states 
are nothing but the surface Andreev bound states. Another case of interest is found 
among the superconductors/superfluids with time-reversal symmetry. In this case 
the edge carries a spin current of the Majorana fermion.SS 
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